Concurrence vs purity: Influence of local channels on Bell states of two qubits 
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We analyze how a maximally entangled state of two-qubits (e.g., the singlet if) a ) is affected by 
action of local channels described by completely positive maps £ . We analyze the concurrence and 
the purity of states qs = £ ®T[%l) a ]. Using the concurrence- us-purity phase diagram we characterize 
local channels £ by their action on the singlet state ip„. We specify a region of the concurrence- vs.- 
purity diagram that is achievable from the singlet state via the action of unital channels. We shown 
that even most general (including non-unital) local channels acting just on a single qubit of the 
original singlet state cannot generate the maximally entangled mixed states (MEMS). We study in 
detail various time evolutions of the original singlet state induced by local Markovian semigroups. 
We show that the decoherence process is represented in the concurrence-us. -purity diagram by a 
line that forms the lower bound of the achievable region for unital maps. On the other hand, the 
depolarization process is represented by a line that forms the upper bound of the region of maps 
induced by unital maps. 



I. INTRODUCTION 

From two well-established properties of the entangle- 
ment, namely, from the fact that i) interactions create 
entanglement and from the fact ii) entanglement cannot 
be shared freely (monogamy [1-4]), we can conclude that 
any non-unitary evolution of a single qubit that is en- 
tangled with another qubit is accompanied with a dete- 
rioration of the original entanglement between these two 
qubits. 

The aim of this paper is to address the question how lo- 
cal actions (channels) affect properties of quantum states 
of bi-partite systems. In particular, we will analyze in de- 
tail how the entanglement and the purity of a two-qubit 
system that has been originally prepared in a maximally 
entangled Bell state depend on the action of a singlc- 
qubit channel, i.e. we assume that one of the qubits of 
the original Bell pair is affected by an environment. 

We have a two-fold task in front of us: First, we will 
analyze how local channels affect the entanglement and 
the purity of the original Bell state. Second, we we study 
how the time evolution (i.e., a one parametric subset St 
of the set of all completely-positive maps) can be rep- 
resented as a one-parametric curve in the concurrence 
vs. purity "phase" diagram. We will focus our atten- 
tion on Markovian evolutions, i.e. those one-parametric 
subsets of channels, for which the semigroup property 
£t£ s — £t+s holds. We will analyze in detail physical 
processes such as decoherence, decay, quantum homog- 
enization, etc. in terms of the concurrence vs. purity 
phase diagram. 

Let us first define those quantities that we shall use 
through the paper. The purity (or equivalently the 
"mixedness") of a state that is described by the density 
operator g will be quantified by the function 



P(g) = Tr[g 2 } , 



(1.1) 



minimum for maximally mixed state, i.e. for the total 
mixture g = 4/, the purity achieves the minimal value 
that is equal to 1/d. 

The entanglement between two quantum systems de- 
scribed by a density operator gAB = Q will be quantified 
by the function called the tangle 



r(e) = 



qkS 2 (ipk) , 



(1.2) 



where tf>k denotes the projection onto a pure state \ipk). 
The minimum in Eq. (1.2) is taken over all pure-state de- 
compositions of the state g while the function S%{iph) = 
2[1 — P(TrsVfe)] is the- so-called linear entropy. The quan- 
tity C(g) — y/r(g) (the square root of the tangle) is 
known in the literature as the concurrence. Wootters [5] 
has derived a simple analytic formula for the concurrence 
of two qubits in a state g 



C(q)=2 max{nj} - \ij 



(1.3) 



which equals to unity for pure states and achieves its 



where fj,j are square roots of eigenvalues of the matrix 
R = g(o y ® <j y )g*((Ty ® o-y) and q* denotes the complex 
conjugation of the original two-qubit density operator g. 
From these definitions it is obvious that the entanglement 
and the purity are closely related quantities and that for a 
given two-qubit state g they cannot take arbitrary values. 

One of the questions one can ask at this point is: 
Which two-qubit states are maximally entangled provid- 
ing that their purity is fixed, and vice versa? This prob- 
lem has been addressed in several earlier papers [6-11]. 
In particular, Ishizaka and Hiroshima [6] have introduced 
the so-called maximally entangled mixed states (MEMS). 
These are the states that for a given value of the purity 
achieve the maximal entanglement. In Ref. [7] a slightly 
more general problem has been solved. The authors have 
shown which unitary transformation has to be applied on 
a given state g in order to maximize the entanglement. In 
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other words: which state maximizes the entanglement for 
a given spectrum of the density operator (i.e., for a given 
value of the purity). After these introductory papers 
have appeared many different aspects of the relation be- 
tween the entanglement and the mixedness have been an- 
alyzed. In addition, the entanglement-based "ordering" 
(parametrization) of the state space of two-qubits, orig- 
inally introduced by Eisert and Plenio in [12], has been 
investigated in detail. It has been shown, that different 
entanglement measures define different ordering of states 
[13]. In fact, this feature is not only characteristic for en- 
tanglement measures, but also for different measures of 
mixedness. This means that the choice of the measures 
affects the final entanglement-purity picture of the state 
space. In Rcf. [8] the analysis of the entanglement-purity 
dependence for various measures has been analyzed. 

In this paper we will study the entanglement-purity 
relation from a perspective of local operations. In what 
follows we will assume that the initial state of two qubits 
is a maximally entangled pure state, i.e. the two qubits 
are prepared in a Bell state. Without the loss of general- 
ity we can consider that the two qubits are prepared in 
the singlet state 



tps = ® I - o x <g> a x - o-y 



i Oy - a z <g> a z ) . (1.4) 



This state is transformed under the action of a completely 
positive trace-preserving linear map £ [14, 15] that de- 
scribes the most general quantum process into the state 



Q£ = £® I[lp s ] . 



(1.5) 



In general, any local action £ (except for unitary opera- 
tions) decreases the purity of the singlet state. Our aim 
is to find how much the entanglement is changed under 
the action of the map £ <g> I. 

In Sections II-IV we will analyze in detail the action 
of unital channels (i.e. those channels that do not affect 
the total mixture, i.e. £[I] = I). The unital channels 
are defined in Sec. II. In Sec. Ill we present a geomet- 
rical representation of the space of all unital maps. In 
Sec. IV we introduce the concurrence- ^s-purity phase di- 
agram and we determine the region that is covered by the 
states qs = £ ®Z[ip s ] that are obtained via the action of 
local unital maps on the singlet state of two qubits. 

The Section V is devoted to investigation of the action 
of non-unital channels. Finally, in Section VI we will dis- 
cuss the time evolution in the concurrence-purity (C-P) 
phase diagram for specific quantum processes. In Con- 
clusion we will summarize the main results and discuss 
some open problems. 



II. UNITAL CHANNELS 

Let us consider firstly the unital channels. Thanks to 
a seminal work of Ruskai et al. [16] the investigation 



(parametrization) of single-qubit channels can be signifi- 
cantly simplified. A single qubit channel £ (not only uni- 
tal ones) can be written as a sequence of two unitary rota- 
tions and one specific completely positive map <&£ which 
belongs to a 6 parametric family of maps. In particu- 
lar, £[g] = U^siVgV^W. Due to the fact that unitary 
operations preserve essentially all interesting properties 
of the original channel £, one can reduce the analysis 
of single-qubit channels into the investigation of proper- 
ties of the channel <&£. In other words, up to a unitary 
equivalence the original 15 parametric set of single-qubit 
channels can be reduced into a 6 parametric set of chan- 
nels $5. This reduction significantly simplifies analysis 
of single-qubit channels. In the Bloch sphere picture the 
general channel £ transforms the Bloch vector f in an 
affine way, i.e. r — > f = Tr + t. The corresponding map 
<&£ acts as follows 



f = Df+ f , 



(2.1) 



where D = diag{Ai, A 2 , A 3 } is a diagonal matrix of sin- 
gular values of the matrix T and f = Rut with Rjj being 
a three-dimensional rotation associated with the unitary 
transformation U . The vector f represents the shift of the 
total mixture, i.e. it is associated with the non-unitality 
of the channel under consideration. 

A set of unital channels form a three-parametric 
family of completely positive (CP) maps of the form $f = 
diag{l, Ai, A2, A3} and the inequalities 



1 + A, - ^ 



A„-A z > 

l-\ x + \ y -\ z > 

1 - A x - X y + \ z > 

l + \ X +\y + \ Z > 0, 



(2.2) 



guarantee the complete positivity of these maps. 

Our task is to evaluate the purity and the concurrence 
of states fig = <i>£ ® I{4> s ] (with ip s = ^ ( J ® J — a x ®a x — 
~y® <J y — °~z®o~z)) as functions of these three-parameters. 
The state Slg takes a simple form il £ = j(I — \ x o- x ®a x — 
)a z ). The corresponding matrix of this 



\y<Jy 



)<Ty—\ Z <J Z 



state reads 



fig = 



/ A D 

B C 

C B 

\D A 



(2.3) 



with A = i(l 



\ z ),B = \{l + X z ),C = -\{\y + \ x ) 
and D = \{\ y - X x ). 

In order to evaluate the purity of the state Q,£ we have 
to find eigenvalues of the matrix (2.3). These eigenvalues 
are given by the expression ki,2 = A ± D and ^3,4 = 
B ±C. Thus, for the purity of the state fig we find the 
expression 



P(n £ ) = Tr[0|] 



(l + A^+A^+A^). 



(2.4) 
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In order to find the concurrence of the state fig we 
have to evaluate the eigenvalues of the matrix 



R = figfTj, ® <T V Q,£<7 



y i £g (Ty 



(X Y 

P Q 

Q P 

V Y X 



(2.5) 



withX = A 2 + D 2 , Y = -2AD, P = C 2 + B 2 , Q = 2CB. 
The square roots of these eigenvalues are fii : 2 = \B ±C\ 
and /U34 = \A ± D\. Since the eigenvalues of fig are 
positive, the square roots of eigenvalues of R and the 
eigenvalues of fig coincide, i.e. the absolute values can 
be removed and the concurrence is given by the formula 



c(fi £ 



X x + Aj, + X z — 1 
An — Aj, — X z — 1 
— X x + X y — X z — 1 
— A K — Aj, + A z — 1 



(2.6) 



III. PARAMETRIZATION OF LOCAL CP 
MAPS: GEOMETRIC PICTURE 

We have derived explicit relations for purity and con- 
currence of a two-qubit density operator fig that is ob- 
tained via the action the unital channel £ on the singlet 
state (1.4). Unfortunately, the concurrence is not so easy 
to deal with. Let us illustrate the whole situation in the 
space of parameters X x ,X y ,X z . It is a well known re- 
sult of the analysis of qubit channels that unital channels 
(specified by A's) form a tetrahedron with unitary Pauli 
operators in its vertices. The same geometrical picture 
holds for states of the form fig. The extremal points 
of this tetrahedron are mutually orthogonal maximally 
entangled states (the Bell basis). The convex combina- 
tions of these states form the tetrahedron, i.e. a classical 
probability simplex. 

The states of the same purity correspond to a sphere 
of the radius proportional to a specific value of the purity 



centered at the point X x = 



A 2 



0. In particular, 



|A| = 4P— 1. There are only four pure states represented 
by the intersection of the tetrahedron with the sphere of 
the radius |A| = 3, i.e. the sphere in which the tetra- 
hedron is embedded. These points are exactly the four 
maximally entangled states forming the Bell basis. 

The equally entangled states specify planes (that form 
a polytope inside the tetrahedron). Because of the dis- 
crete symmetry represented by four unitary transforma- 
tions (rotations) described by the operators /, a x ,<7 y ,a z , 
the analysis of possible values of the concurrence and the 
purity in terms of A can be focused onto only two cases: 
i) all Xj are positive, or ii) all Xj are negative. The sym- 
metry relating these two options is the space inversion 
(A — ► — A), which is not physical (i.e., this inversion can- 
not be realized by a CP map). In fact, the tetrahedron 
does not possess such symmetry. The tetrahedron is sym- 
metric under rotations by an angle <j) — 7r/2 along each 







^1,1,1) 











"(-1,1,-1) 



(1,-1,-1) 

FIG. 1: (Color online) This figure presents the space of 
unital channels <E>g parameterized by A = (Ai,A2,As) (with 
— 1 < Aj < 1). Up to unitary rotations the set of states Q,£ 
for that are obtained by the action of unital channels £ on 
the singlet state (1.4) form a tetrahedron with vertices cor- 
responding to maximally entangled states that form the Bell 
basis. The octahedron inside the tetrahedron represents a set 
of all separable states. The remaining four tetrahedrons in 
each corner of the original tetrahedron represent sets of en- 
tangled states. The closer the point is to the vertex the more 
entangled the corresponding state is. The planes parallel to 
faces of the separable region contain states with the same de- 
gree of entanglement ( "iso-concurrence" planes). States with 
the same purity form spheres ( "iso-purity" spheres) centered 
at the center of the cube, i.e. at the point A = (0, 0, 0). 



axis (cr matrices). Consequently, the analysis of the max- 
imum of the concurrence reduces to an investigation of 
two cases (see Fig. 1): Xj > (for all j) and Xj < (for 
allj). 

Let us start with the case of positive values of As. 
In this case the maximum of the concurrence (2.6) is 
achieved for C(fig) = \{X X + X y + X z — 1). Providing 
that this number is larger than zero, then states of the 
same concurrence C form the plane X x + X y + X z — d with 
d = 2C + 1 . These "iso-concurrence" planes intersecting 
the tetrahedron (its positive octant) are given by normal 
vectors n = (1, 1, 1), i.e. ft ■ X = 2C + 1. The plane 
with C = (d = 1) form a boundary (face) of separable 
states in this part of the tetrahedron. Due to the sym- 
metry mentioned above the same picture holds for other 
four vertices of the tetrahedron, i.e. tetrahedrons "un- 
der" four maximally entangled states (that correspond 
to vertices of the tetrahedron) . 

The negative region of the tetrahedron (Aj < 0) con- 
tains no entangled states, i.e. it consists of only separa- 
ble states. One can prove this by analyzing all possibil- 
ities, or also by exploiting the geometrical picture. The 
set of allowable (negative) A is bounded by the plane 
1 + X x + X y + X z =0 that potentially contains entangled 
states in this "negative region". If not, then due to a 
convexity of separable states the whole "negative region" 
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is separable. Using the expressions for X x ,X y ,X z deter- 
mined by the equation of the plane and by calculating the 
concurrence we obtain C = \ max{— 2, 2A X , 2A a , 2A 2 , 0}. 
Because of Xj < we obtain that the concurrence always 
equals to zero. 

So far, we have shown that entangled states belong to 
regions close to the vertices of the tetrahedron of all pos- 
sible states f2g. Separable states form the octahedron 
embedded in the tetrahedron of all states (see Fig. 1). 
Our next task is to evaluate the concurrence and the 
purity for all entangled states. In particular, we are in- 
terested in the shape of the region formed by the states 
fig in the C-P phase diagram. 



IV. CONCURRENCE VS PURITY DIAGRAM 

One possibility how to characterize bi-partite states is 
to specify their purity and the value of their entanglement 
[6]. Such classification contains highly non-trivial infor- 
mation about the state itself. One can choose different 
measures for both the purity and the entanglement. The 
resulting characterization strongly depends on the partic- 
ular choice of measures of entanglement and the purity 
[7, 8] . In what follows we will use two most common mea- 
sures: the mixedness for the purity and the concurrence 
for the entanglement. 

As we have already mentioned the boundaries of the 
C-P phase diagram has been analyzed and the results 
are known. The states maximizing the concurrence for 
the given value of purity are known [6] as maximally en- 
tangled mixed states (MEMS). This concept generalizes 
the notion of Bell states, i.e. maximally entangled pure 
states. MEMS have the following form (up to local uni- 
tary transformations): 

&MEMS = p\ 4 > +) (.4*+ 1 + (l-p)|01)<01| (4.1) 
for p G [2/3, 1] and 

Pmems = p|0+)(^+| + i|Ol)(Ol| (4.2) 
+ (| -p/2)(|00><00| + |11)<H|) 

for p G [0,2/3] and [<£+) = ^(|00) + |11)). These states 
specify the region in the C-P phase diagram which is 
physical, i.e. any point in this part of the C-P phase 
diagram corresponds to a state of quantum-mechanical 
system (see Fig. 2). Our aim is to analyze this picture 
for the states of the form fig, i.e. those states that are 
obtained by the action of local unital operations on the 
maximally entangled state ip s . In particular, our task is 
to find the maximum/minimum value of the concurrence 
for the given value of the purity. 

As we have already shown states of equal purity cor- 
respond to the sphere (parameterized by X x + X y + X 2 Z = 
4P — 1) in the tetrahedron of all states fig (see Fig. 1). 




0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

purity 



FIG. 2: (Color online) The concurrence vs. purity phase di- 
agram of two qubits. The region of states tts — £ <8> I[ips] 
(region colored in red in the C-P phase diagram) that are 
obtained from the maximally entangled two-qubit states by 
the action of unital channels £ is bounded by two lines cor- 
responding to Cw and Cd for a given value of the purity. 
The upper bound of the region of physically realizable states 
(shaded sector of the C-P diagram) corresponds to maximally 
entangled mixed states (MEMS). We note that the MEMS 
states cannot be achieved from Bell state by applying local 
channels on a single sub-system only. 



States with the same amount of entanglement deter- 
mine a plane X x + X y + X z = 2C + 1. We have argued 
that it is sufficient to consider only the positive region 
(Xj > 0) of the tetrahedron. Given the iso-purity sphere 
the iso-entanglement plane with the maximal value of en- 
tanglement is the one that "intersects" the sphere only 
in a single point. One can exploit the geometric pic- 
ture to see that this point (state) fulfills the condition 
X x = X y = X z = X. After we insert this condition into 
the equations for the purity and the concurrence we ob- 
tain 

o\2 zip — i 1 1 , 

3A = 2C + 1 } C ™ = 2 (^4P^i) - 1) (4.3) 

These states form the upper bound of the available region 
in the C-P phase diagram (see Fig. 2). In particular, 
this bound contains (is formed by) Werner states, i.e. 
gw = qips + (1 - q)\l and C max = C W - 

The next step is to analyze the minimum of the concur- 
rence for a given value of purity. In general this minimum 
is trivial, because there always exist separable states for 
a given value of purity. However, in our case we inves- 
tigate states that are generated by local unital channels 
from the maximally entangled Bell states. In this case a 
nontrivial lower bound exists. Exploiting the geometry of 
the states fig we can conclude that this minimum is zero 
for all purity spheres for which the intersection with the 
tetrahedron contains the plane Ai + A 2 + A3 = 1 (C = 0), 
i.e. the boundary of the separable states. The states 
of the maximal purity belonging to this plane (i.e., the 
maximally mixed separable states) are the points with 
two same components and the remaining one equals to 
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unity. Without the loss of generality let us consider the 
case X z = 1 and \ x = \ y = A. The equation of the 
plane implies that A = and consequently, the purity 
P = |(1 + |A| 2 ) = \. That is, for the states with the 
purity larger than 1/2 there are no separable states 

The question on the minimum value of entanglement 
for a given value of the purity is equivalent to the question 
about the maximum of the purity for a given value of 
the concurrence. Using the same arguments as before we 
find that the state of the maximal purity satisfy the same 
conditions, i.e. \ z = 1 and \ x = \ y = A. Consequently, 
the equation of plane \ x +X y -\-\ z = 1+2C implies A = C. 
As a result we obtain that P = |(1 + C 2 ). Inverting 
this formula we obtain the functional dependence of the 
concurrence as a function of the purity that specifies the 
lower bound of the allowable region in C-P phase diagram 
for states fig 



C min = V2P~T 



(4.4) 



In conclusion, states induced by local unital channels 
from the initial maximally entangled state are repre- 
sented in the concurrence vs purity diagram by points 
in the region that is bounded from above by the "line" 
C max = Cw and from the below by the line C m ; n (see 
the red region in Fig. 2). The whole region of physically 
relevant states is determined by the line Cmems corre- 
sponding to the MEMS. Unital channels do not allow 
us to achieve states that are within the lines Cw and 
Cmems • Simultaneously, we can conclude that the states 
fig with the purity larger than 1 /2 remain entangled (this 
property is given by the lower bound C m ; n = Co). 



V. NON-UNITAL CHANNELS 

In the previous section we have shown that unital chan- 
nels (up to unitary rotations) are characterized by three 
parameters. The non-unital channels arc characterized 
by six parameters. The triple f = (t x ,t Vi t z ) describes 
how the total mixture is transformed under the action of 
non-unital channels. In particular, under the action of a 
general non-unital map £, the maximally entangled state 
(e.g., the singlet) is transformed into the state 



t 4 



(7 + f ■ a) <g> I - A • (a <g> a) 



(5.1) 



where A • (a <g> a) = \ x cr x <g> cr x + \ v o y ® u y + A z cr z <g> a z . 
In other words 



fA+r F D \ 

F* B-t C 

C B+t F 

V D F* A-t ) 



(5.2) 



where A, B, C, D are defined as before and F = (t x — 
ir y )/A, t = t z /A. First we address the question whether 
by using non-unital channels the upper bound on con- 
currence can be increased, in particular, whether fig can 



be a MEMS state (i.e. the state on the Cmems line). 
In order to answer this question we use the identity 
Tr^Qf = TiAtps = \l which holds for any local action 
£ ® J. For MEMS states we find that Tta£mems ^ \l 
and Tr^gMEMS ^ \l as well, i.e. Slg cannot be a MEMS 
state. Therefore, we conclude that the MEMS states can- 
not be achieved from the maximally entangled state via 
local operations, i.e. MEMS are not the states of the 
form n £ = $ £ (g)I[ip s }. 

However, the action of non-unital channels on maxi- 
mally entangled states can be different then the action 
of unital channels. That is, the states fls = $g ®T[i(> s ] 
that are obtained by the action of non-unital channels 
might lie in the region of the concurrence vs purity 
phase diagram that are bounded by the lines Cmems 
and Cw- The purity of states fig can be calculated by 
finding the trace of A direct calculation gives that 
P = i(l + |A| 2 + |r| 2 ). Thus in the picture of A param- 
eters (for a fixed vector f) the states of the same purity 
belong to the sphere |A| 2 = AP - 1 - |f| 2 . Note that 
for the fixed f the set of all possible A do not form the 
tetrahedron anymore. 

Let us consider a specific case r = (0, 0, r 2 ), i.e. F = 0. 
In this case the eigenvalues of fig read: 



Mi = 


J(l-A,+ 




-A,) 2 + r|) 


m = 






-A,) 2 + r|) 


M3 = 






+ A,) 2 + r 2 ) 


fn = 


J(l + A,- 




+ A,) 2 + r 2 ) 



The positivity of these eigenvalues determines the set of 
all allowed values of parameters A, r z . Fixing t z we obtain 
four surfaces (setting fij = 0) in the three-dimensional 
space of parameters A that form boundaries of the set 
of all possible states fig. The identities fij = can be 
rewritten into the equations 



\ z = 1 



+ A, 



) 2 +r 2 



A, 



-1 



(5.3) 



that completely specify the shape of the set in the space 
of parameters Ai, A 2 , A 3 . For non-unital channels the cor- 
ners of tetrahedron are smoothed ( "rounded" ) depending 
on the value of the shift f. 

To compute the concurrence C of state fig = $£<g>Z[i/> s ] 
we have to find eigenvalues of the matrix 



R = £l£(<Jy ® a y )fl £ (<Ty ® Uy) 



( * 




\5- 



s + 

(3 7 - 

7+ P 

0a 



where a 



A 2 



2C(B ± r), and 8± 



D 2 , (3 = B 2 - 
2D{A ± t) (t 



+ C 2 , 7± = 
z/A). Square 
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roots of the eigenvalues of the matrix R read 

{VA 2 - r 2 ± \D\, VB 2 - t 2 ± \C\} , (5.4) 

and for the concurrence we obtain the expression 

C = max{0, 2{\D\ - \J B 2 - r 2 ), 2(|C| - ^A 2 - r 2 )|5.5) 

Using the parameters A, t z the concurrence can be rewrit- 
ten in the form 

C=imax{ \\-\ x \- y fl + A,) 2 -r 2 2 , 1 . (5.6) 
[ |A, + A,|- V(l-A z ) 2 -r| J 

As we will show in the next Section the states fig = 
$£ (g> Iff/'s] that are generated by non-unital channels 
can lie above the line C ma x in the concurrence vs purity 
diagram. In order to have more physical insight into 
the action of the local channels on maximally entangled 
states let us consider one-parametric set of local maps 
that correspond to specific time evolutions of two-qubit 
systems. 



Due to the evolution of one qubit the singlet is trans- 
formed into a state f2 t with the purity 

^|(l + M 2 )^(l + e - 24/T ) (6.3) 

and with the concurrence 

C t = l(Ai(t) + A 2 (t) + A 3 (i) - 1) = e- t/T • (6.4) 

Comparing these two equations we find that the purity 
and the concurrence of the state f2 t = &£®T{ip s ] induced 
by decoherence acting on one qubit are related as 

Pt = l(l + C?), (6.5) 

or equivalently 

C t = V^Pt - 1 • (6.6) 

Since P t € [1/2, 1] we have obtained that the process of 
single-qubit decoherence in a C-P diagram corresponds to 
the lower bound of the allowed region for unital channels, 

i.e. Ct ■= Cd = C m in- 



VI. EVOLUTION IN C-P DIAGRAM 



B. Depolarization 



In this section we will analyze how the evolution of 
a maximally entangled state under the action of a lo- 
cal channel is reflected in the C-P diagram. The case 
of unitary evolution is trivial: Under the action of lo- 
cal unitary transformations neither the concurrence nor 
the purity are changed. Therefore, the state fit = 
(Ut ® I)ip s (U-t <8> I) is still a maximally entangled pure 
state. In what follows we will analyze several models of 
non-unitary dynamics. We will focus our attention on 
Markovian semigroup dynamics. In particular, we will 
consider processes of the decoherence, the depolarization 
and the homogenization. 



A. Decoherence 

The decoherence of a qubit is induced by the master 
equation [17] q = i[H,g] + (T/2)[H,[H, g]]. The solu- 
tion of this equation form a semigroup of unital quantum 
channels 



St = 



(I \ 

e~'/ T cos ut e~*/ T sinu>i 

— e~*/ T sinw£ e~*/ T coswt 

\ 1 J 



(6.1) 



The singular values of these channels are 



Ai(t) =X 2 (t) = e 
A 3 (t) = 1 . 



-t/T . 



(6.2) 



The process of a single-qubit depolarization in a spe- 
cific basis is represented by the semigroup [14] 



( 1 

e-*/ T 


\ o o 



-t/T 









-t/T 



\ 



(6.7) 



In this case the states f2i are Werner states, i.e. this 
type of dynamics in C-P diagram is represented by the 
line Cw for unital channels. In particular, 



ft, 



-t/T 



(6. 



Thus, we have found two processes that saturates the up- 
per and lower bound of the region in C-P diagram that 
is allowed for unital channels. In particular, the decoher- 
ence saturates the lower bound while the depolarization 
process defines the upper bound. 



C. Homogenization 

The process of quantum homogenization is described 
by the semigroup of non-unital channels [18, 19] 



-t/T-, - 





e-^R^t 








-t/T 



(6.9) 
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where 



cos u>t sin ult 
— sinwi cos ujt 



(6.10) 



is a rotation matrix. This process describes an evolu- 
tion that transforms the whole Bloch sphere into a single 
point, i.e. a generalization of an exponential decay. That 
is, quantum homogenization is described by a contractive 
map with the fixed point that is the stationary state of 
the dynamics. The parameters in the description of the 
map (6.9) have the following meaning: w is the purity 
of the final state, T\ is the decay time, Ti is the deco- 
herence time, and lo describes the unitary part of the 
evolution. The singular values are similar to those in the 
decoherence, i.e. 



A 1 (t)=A 2 (i)=e-*/^ 
Aa(t) = e-*/ T i . 



(6.11) 



The homogenization belongs to a class of non-unital 
channels that we have analyzed in previous section. 
Therefore we can easily find expressions for the purity 



1 



1 + 2e" 2t / T2 + e" 2 */ Tl + u- 2 (l - e- { / Tl ) 2 j6.12) 



and the concurrence 

C t = max{0,e-* /T2 - i(l - e~ t/Tl )y/l-w 2 } .(6.13) 

The resulting lines for all considered evolutions (decoher- 
ence, depolarization and homogenization for the values 
Ti/T 2 = 1/2 and w = 1) are depicted in Fig. 3. 

In a special case, when w — (final state is the total 
mixture), the homogenization process is unital. In this 
case the purity and the concurrence are 

P t = |(1 + e- 2t / Tl + 2e~ 2t ^) (6.14) 
C t = e~*/ T2 + i(e~*/ Tl - 1) (6.15) 

In Fig. 5 we can see the corresponding line for different 
fraction of decay and decoherence times, i.e. T\ and T 2 , 
respectively. Interesting point is when these two rates co- 
incides (Ti = T2) when the homogenization saturates the 
Werner states line (i.e. CV). In the limit of X1/T2 — * 00 
the homogenization approaches to the decoherence line 
(i.e. C D ). 

For w = 1 the homogenization describes the exponen- 
tial decay to a pure state. In such case 



Pt = 1(1 



e -2t/Ti + e -2t/T 2 

C t = e"^ 



(6.16) 
(6.17) 



In this case one can express the purity as a function of 
concurrence, i.e. P = 1(1 + C 2 + C 2T ^ - CF^). In 
the special case T 2 = 2Tj (i.e. decoherence time is twice 
as fast as decay time) the homogenization follows the line 



. 8 
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C = \/2P-l 



(6.18) 



FIG. 3: (Color online) Parametric plots of time evolution of 
the concurrence and the purity in processes of the decoher- 
ence, the homogenization and the depolarization. The lines 
are parameterized in such way that the initial moment of the 
time evolution t = corresponds to the point C = P = 1 
(upper right corner of the C-P diagram). The decoherence 
line (2) represents the lower bound (Cd) of the unital region, 
the depolarization line (1) forms the upper bound (Cw) of 
the unital region. The homogenization (a non-unital process) 
is characterized by the line (3) that is outside the unital re- 
gion. In the present case the homogenization describes an 
exponential decay, i.e. the fixed point of the evolution is the 
state jO). 



This line goes above the region determined by unital 
channels (see Fig. 3). 

In Figs. 4-6 we analyze various regimes of the homoge- 
nization process: We consider processes with a fixed value 
of the ratio J1/T2 = 1/2, but we change fixed points of 
the evolution (Fig. 4). In Fig. 5 we consider the homog- 
enization process with the fixed point being equal to the 
total mixture (i.e. in this case the homogenization is a 
unital process) and we consider different values of the 
ratio Ti/T 2 . In Fig. 6 we consider the homogenization 
process with the fixed point equal to a pure state (i.e. in 
this case the homogenization is a non-unital process) and 
we consider different values of the ratio T\jT%. 



VII. CONCLUSION 

In this paper we have investigated how a local trans- 
formation (quantum channel described by a CP map) 
of a sub-system affects the entanglement and the global 
purity of the initial singlet state (or, equivalently, ar- 
bitrary maximally entangled states) of two qubits. We 
have analyzed the states gg — £ ®I[i/j s ]. We have used 
the concurrence and the purity of these states to classify 
local channels £. In particular, using the concurrence- 
US. -purity diagrams we have specified the region induced 
by local unital channels. This region does not cover the 
whole set of physically realizable states (see Fig. 2). Lo- 
cal unital channels induce states that are represented by 
points in the region of the C-P diagram that is bounded 
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purity 



FIG. 4: (Color online) We present evolutions of the concur- 
rence and the purity for the homogenization process in the 
case with the ratio of characteristic times Ti and Tb taking 
the constant value T1/T2 = 1/2 while the fixed point of the 
evolution varies from total mixture (unital process) to a pure 
state (non- unital process) . In particular, when the fixed point 
of the homogenization is a total mixture then the homogeniza- 
tion dynamics in the C-P diagram is described by the line (1). 
This line lies below the line Cw with the final point C = 
having the smallest value of the purity. On the other hand, 
when the fixed point of the evolution is a pure state, then the 
corresponding C-P line (2) ends in the point C = with the 
largest value of the purity P = 1/2. 




0.3 0-4 - a n.fi 0.7 0.8 0.9 1 



purity 

FIG. 5: (Color online) The evolution of concurrence and pu- 
rity for the homogenization with the fixed point equal to the 
total mixture. In this case the homogenization is a unital 
process. We vary the ratio T1/T2 in the interval [1/2, 00]. 
For Ti = T2 the homogenization exactly covers the line (2) of 
Werner states (Cw)- In the limit of T1/T2 — > 00 this evolution 
coincides with the line of decoherence (Cd)- 



from below by the line Cd and from the above the line 
Cw . We have shown that even the most general (includ- 
ing non- unital) local channels acting just on a single qubit 
of the original singlet state cannot generate MEMS. This 
means that the upper bound specified by the Cmems 
line cannot be achieved by the action of the local chan- 
nel of the form gg — £ ®T[ip s \. Specific achievable upper 
bound for non-unital maps is to be determined. It def- 
initely is above the line Cw (except the values of the 
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purity 



FIG. 6: (Color online) The parametric plot of the concurrence 
and the purity for the homogenization with the fixed point 
that corresponds to a pure state. That is we consider an 
exponential decay into a pure state. We vary the T1/T2 in the 
interval [1/2, 00]. In this case, the evolution ends in the point 
C — and P = 1/2 irrespective of the particular value of the 
ratio T1/T2. When T1/T2 = 1/2 the evolution is represented 
by the line which is "deep" in the non-unital region of the C-P 
diagram. On the contrary, for Ti/Tb — * 00 the corresponding 
line coincides with the bound Cd (the decoherence line). 



concurrence C — 1 and C — when it coincides with the 
bound on unital maps) . From our analysis it follows that 
for unital maps the lower bound of the achievable region 
is determined by the line Cd- We conjecture, that this 
is also a general lower bound for non-unital maps. This 
conjecture is supported by our numerical analysis, but 
we have no rigorous proof yet. 

From our previous analysis an interesting observation 
follows. Specifically, if the state gg — £ <£> Z[ip s ] has a 
purity larger than 1/2 then the concurrence has to be 
non-zero, correspondingly, the state is entangled. On the 
other hand, if the purity is less than 1/3 then the state 
is separable. 

In our paper we have analyzed neither the action of 
the bi-local channels (£1 <g> £2) nor the action of the non- 
local maps £\2- It is clear, that using non-local maps any 
point of the C-P diagram below the MEMS line Cmems 
can be achieved, i.e. the state gs = Swtys] can have ar- 
bitrary value of the concurrence and the purity that is in 
the region specified by the bound Cmems- On the other 
hand, bi-local operations generate states gs = £1 (S^ft^s] 
that are represented by the points in a region of the C-P 
diagram that is restricted from above and from below. 
Specific boundaries are not known. From particular ex- 
amples one can conclude that if the two local transforma- 
tions are non-unital maps representing the exponential 
decay then one can achieve states with zero concurrence 
but maximal purity. We will analyze these bounds else- 
where. 

We have studied in detail time evolutions described 
by Markovian semigroups. In particular, we have shown 
that the decoherence process is represented by a line that 
forms the lower bound Cd of the achievable region for 
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unital maps. In our parametrization the time t = is rep- 
resented by the point C = 1 and P = 1, while the point 
of the entanglement destruction (C — 0) is achieved at 
some "entanglement-breaking" time t sep that is infinite - 
see discussion below. 

The depolarization process saturates the upper bound 
(Cw) of the region of maps induced by unital maps. 
Here the entanglement-breaking time is finite and equals 
t S ep — Tin 3. For dynamics that are described by 
non-unital maps (e.g. homogenization processes) the 
entanglement-breaking times can be both finite as well 
as infinite (as in the case of the exponential decay). 

We have paid attention only to evolutions governed by 
semigroups, i.e. by Markovian processes. A general rule 
is that for this type of time evolutions the associated C-P 



line must be non-increasing, because local action cannot 
create entanglement. This property prevents from loops 
in C-P diagram, but still allows that the purity increases 
while the concurrence is decreasing (see Fig. 6). For non- 
Mar kovian evolutions it is possible to observe even loops, 
however, the question of more general dynamics is out of 
the scope of this paper and will be discussed elsewhere. 
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